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Abstract 

Recently the quantum nature in the energy transport in solar cell and light-harvesting com¬ 
plexes have attracted much attention, as being triggered by the experimental observations. 

We model the light-harvesting complex (i.e., PEB 50 dimer) as a quantum heat engine (QHE) 
and study the effect of the undamped intra-molecule vibrational modes on the coherent en¬ 
ergy transfer process and quantum transport. We find that the exciton-vibration interaction 
has non-trivial contribution to the promotion of quantum yield as well as transport properties 
of the quantum heat engine at steady state, by enhancing the quantum coherence quantified 
by entanglement entropy. The perfect quantum yield over 90% has been obtained, with the 
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exciton-vibration coupling. We attribute these improvements to the renormalization of the 
electronic couplings effectively induced by exciton-vibration interaction and the subsequent 
delocalization of excitons. Finally we demonstrate that the thermal relaxation and dephasing 
can help the excitation energy transfer in PEB 50 dimer. 


Introduction 

Recently the wide-spread interest in uncovering the quantum phenomena in the solar cell and pho¬ 
tosynthetic process has been triggered by the experimental investigations of the electronic dynam¬ 
ics in light-harvesting and Fenna-Matthews-Olson (FMO) complexes.-- The transport of excita¬ 
tion energy absorbed by antenna towards the reaction center (RC) occurs with a perfect efficiency 
over 90%,-» 6 which was shown to be strongly correlated to the long-lived quantum coherence be¬ 
tween different molecules.- 7 ’ 8 

Numerous research made it clear that the exciton energy transport in photosynthetic organism 
critically depends on the exciton-phonon interaction,- - - 2 besides the electronic coupling between 
molecules. The exciton-phonon interactions have two types: high-frequency modes from the nu¬ 
clei vibrations in molecules, and low-frequency modes induced by environmental fluctuations. 
Actually these two kinds of exciton-phonon interactions in these complexes are associated with 
low-energy fluctuations of protein immersed in the solvent and intra-molecular vibrations, 13114 
respectively. Since the intra-molecular vibrations are undamped, it is evidently shown to have sig¬ 
nificant influence on the coherent energy transfer when energy quanta of vibrational modes is in 
resonance with the energy splitting of excitons.- 15- *^ 

In this work, we study the effect of intra-molecule vibrations on energy transfer processes in 
the pairs of chromophores, which describes several light-harvesting antennae in the nature. Two 
important examples are the central PEBsoc — PEB^/ dimer in the cryptophyte antennae PE545 
(Phycoerythrin545) and the Chicot ~ Chl fl 602 pair in the light-harvesting complex II (LHCII). 
The natural photosynthetic organism functions in the presence of both solar radiation and low- 
energy fluctuations of protein. In our model, there are three different types of energies involved, 
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namely solar radiation, high-frequency modes from the nuclei vibrations and low-frequency modes 
from the environmental fluctuations. According to the quantum thermodynamics, QHE converts 
hot thermal radiation into low-entropy useful work. 18-20 We will include the dynamics of dis¬ 
crete vibrational modes between intra-molecules, together with the dynamics of system, due to the 
comparable relaxation time. 8,21-22 We will find that the interaction with vibrational modes has a 
non-trivial contribution to the enhancement of the nonequilibriumness flux 23 - 24 trapped by RC and 
energy transfer efficiency (ETE), by improving the coherence. 




Figure 1: (Color online) (top) Illustration of cryptophyte antennae phycoerythrin 545 (PE545); 
(bottom) Schematic of our model for PEBsoc-PEBsod dimer in PE545. Standard parameters are 
taken from Ref.: 29 - 30 ft) = 800 cm -1 being quasi-resonant with the intra-molecular mode of fre¬ 
quency, £i — £2 — 1042 cm -1 , A = 92 cm” 1 and 7? = 300K. 


Polaron transform and quantum master equation 

We consider a prototype dimer where each chromophore has an excited state described by exciton 
with energy £, strongly coupled to a quantized vibrational mode of frequency ft), with the identical 
coupling strength A for i = 1,2. The Hamiltonian of dimer and vibrational modes read 

Hex = -E'o|0)(0| +£icjci +£ 2 4 c 2 +A(c{c 2 + 4 c i) 
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and H v ib = h(o{b\b\ -\-b\b 2 ). The interaction term is 


H- e x-vib — V2Uuoc[ci(bi + b\) + V2Uuoclc 2 {b 2 + b\) (2) 

In what follows we will reduce the degrees of the exciton-vibration dynamics, by introducing the 
correlated (+) and anti-correlated (-) vibrational coordinates: b± — (b \ ±b 2 ) /y/l, which rewrites 
the Hamiltonian as 


Hex-vib = kh(D (cjci — c|c2^ (b-+bl_^j +?ih0)(c\ci+clc2^ (b + + b\^j (3) 

where the 1 st term in H ex - V ib describes the vibration mode with phase difference of K while the 
2 nd term describes the center-of-mass motion of vibrational modes. Thus the anti-correlated vibra¬ 
tional mode is excited by optical phonons.— Next we will employ the polaron transformation— 
to eliminate the center-of-mass motion of intra-molecule vibrations and uncover the effect of the 
exciton-vibration interaction. The generating function can be written as S — A [(cjci — c^c^ibl — 
bJ) + (cjci +c\c 2 ){b' + — b + )\, which is obviously anti-Hermitian, namely, S’*' = —S, to ensure the 
unitarity of the transform. The new Hamiltonian in the polaron frame is obtained by evaluating 
e s He~ s , and it takes the form of 


H = {e 1 — 2A 2 /z<u) c\c 1 + (£2 — 2k 2 hoo) 


(4) 


+ A 


e 2 W-- b -)clc 2 + e- 2 W-- b -)cl Cl 




vib 


where H vi b = hco(b f x _b + +b^_b^). From Eq.(??) it is clear that the center-of-mass motion is in¬ 
dependent of the dynamics of exciton so that we can ignore its effect. The exciton-vibration in¬ 
teraction causes a renormalization of the electronic coupling, which as will be shown later, will 
critically affect the quantum transport. Notice that the one-particle approixmation for exciton is 
not used here and after. 

The Hamiltonian in polaron frame provides a picture that the transition of exciton in the dimer 
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is assisted by exchanging energy with the vibrational modes. Physically we can only consider the 
single-quanta process of vibrations, owing to: (/) the low probability of multi-excitation and (if) the 
energy scale of the exciton-vibration interaction being in quasi-resonance with the gap between the 
adjacent levels of vibrational modes. Hence the dynamics of intra-molecule vibrations can be re¬ 
stricted into the space spanned by { \m ), \m+ 1)}, as shown in details in supplementary information 
(SI). Moreover the occupation m here is the average of particle number in its eigenmodes rather 
than the thermal occupation of bosons. The total Hilbert space is M' = 00^eiec it,. 

In light-harvesting complexes, i.e., PEB 50 dimer and LHCII pair (or FMO complex), the ex- 
citons need to interact with radiation environment from solar as well as low-energy fluctuations 
of protein immersed in the solvent, in order to realize the energy transfer to RC. In addition, we 
need a connection of one site in the dimer (labelled by 2 in our notation) to RC, described by a 
trapping rate F, in order to generate the output work. The interactions to radiation and protein 
reservoirs read (notice that this is not influenced by the polaron transform introduced before, since 
the polaron transform only operates in the subspace of the system Hamiltonian) 

2 2 

flint = £ £#k p Ui + ct) (a kp + al k p ) + £ £ Yq,sficjci (bqs + bl q s ] (5) 

/=lk,p i=l q,i 

where a kp and b qs are the bosonic operators for radiation and low-energy fluctuation reservoirs, 
respectively, p and s denote the polarizations of the modes of radiation and low-energy fluctuation 
reservoirs, respectively. Thus the total Hamiltonian of the system and environments reads H — H + 
//bath T H'mu where //bath = Lk.p a li P a ^p + Lq.s b' qs b qs . Based on the perturbation theory, the whole 
solution to density operator can be written as Psr = PsO) 0' Pr( 0) -t-pg(f) with the traceless term 
in higher orders of coupling. Since the time scale associated with the environmental correlations is 
much smaller than the time scale of the system over which the state varies appreciably, the quantum 
master equation (QME) for reduced density matrix of systems can be derived under the so-called 
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Markoff approximation 


dp 

dt 


7 7 

j[P,fi} + -^2 L L r «^c(l«)(^|p|c)(^| -Sad\c){b\p) 

" n a,b= 0 c,d =0 V 7 

+ i E E I ^ (|a> (% |c><c| - S flC |c><%) + h.c. + r^(p) 

" fl= 0 fc= 0 c =0 V 7 


( 6 ) 


where the decay rate T adbc /h 2 = v dc Lu Lm, n =o r S v mn U h] U ak U n iUj h and the dephasing rate R abc /h 2 = 

Y]=oLk.l Ykt PiP c UikUjjU'l a U b i. The superoperator describing the trapping by RC is % r {p) = Ln=m (|0,0, n) (0,1, n\f 
|0,l,n)(0,l,n|p|0, l,n)(0,1,n|).r[* =/dq/(2^) 2 g2(„ q + 1)5(0^ — Vq ) fork>l,or Jdq/(2K) 2 g\n q S(a^ d + 
v q ) for k < l. For the environment of protein in the solvent, we use the Debye spectral density: 

J(co) — ^ o^+co 1 w ^ crc Er is so-called reorganization energy, and y^f = J(\co'ki\)n(\(o'ki\) for 
k > I, J{\(o'ki\) [n{\(o\i\) +1] for k < l, or for k — l. U is the unitary matrix that di¬ 

agonalizes the Hamiltonian matrix H. In the Liouville space, the QME can be formulated as 
two-component form 


d 



dt 


\P C J 


p pc 

H 

\^S$cp ■ // $c J 

\Pc) 


(7) 


Here, p p and p c represent the population and coherence components of density matrix, respec¬ 
tively. To obtain the population dynamics, we need to project the QME into population space, 
by eliminating the coherence components using Laplace transform.— Based on the view point of 
QHE, the whole system should work at steady state, which is what we are interested in this Letter. 
Thus in the long-time limit, the QME at steady state reads 


( J( p - Jt pc J( c 1 .M cp ) p S p S = 0 ( 8 ) 

where we introduce stf = — .Mp C .M ~ 1 .M cp and can define the transfer matrix as T mn = stf n ,„p£. 

Thus the net-nonequilibrium-flux is c mn = £# nm Pm — min (^nmpm^mnpn), which can be decom- 
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posed into the superposition of several closed loops . 24,27 Actually c mn quantifies the detailed- 
balance-breaking and time-irreversibility. In reality, people are principally able to observe the flux 
trapped by RC, which is defined as JPtrap — + P 5 )• This is consistent with the definition of 

c mn since the transition induced by RC is unidirectional. As we will see later, ^ trap will play an 
significant role in the discussion of transport properties of this QHE. 


Coherent energy transfer and quantum transport 


Energy transfer efficiency and flux trapped by RC 

Now we are able to discuss the coherent energy transfer in the dimer, after the absorption of photons 
from the Sun. In the natural light-harvesting complexes, the reorganization energy is Er = 34cm , 
the cut-off frequency is (O c i = 50 fs -1 , the temperature of the radiation reservoir is 1\ — 5780K— 
and the trapping rate by RC is set to be T — lps“ 1 . In PEB 50 dimer the dipole moments in Eq.(??) 
are f\ = 1.0, — —0.9. 29 First we introduce the energy transfer efficiency (ETE) 10, 11 




trap 


cfitrap T K^g 


(9) 


where the decay rate is Kj e — L), -0 L v=2 ^ vPv — r(p^ + pf), which quantifies the probability 
going back to ground state per unit time. Figj2]X) and Ob) illustrates the effect of coupling strength 
of exciton to intra-molecule vibrations on the trapping flux by RC as well as ETE. As we can see, 
both of the flux and ETE show a sharp increase as the coupling strength becomes large, besides a 
small decrease at beginning. The reason for such decrease is the corresponding small increase of 
population on site 1 as shown in Figllfc), which indicates that the exciton becomes more localized 
at very weak coupling to molecular vibrations. Significantly, the large coupling to vibrations leads 
to the optimization of ETE of the system, namely, over 90%. Physically, we can understand it as 
follows: the exchange of the energy between excitons and vibrational modes leads to the renor¬ 
malization of the electronic coupling, which effectively amplifies the magnitude of the electronic 
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couplings by a factor of y/m+ 1 as shown in the Hamiltonian. Consequently, the exciton trans¬ 
port is accelerated. On the other hand, Figj2}a) also indicates that the nonequilibriumness of the 







Figure 2: (Color online) (a) Flux trapped by RC, (b) energy transfer efficiency, (c) the population 
on site 1 (small figure) & the square root of ratio of the populations of two chromophores in the 
PE545 dimer (large figure), (d) the entanglement entropy (which also quantifies the coherence 
effect), (e) work generated by QHE to RC and (f) entropy production (EPR) vary as a function of 
coupling strength of exciton to intra-molecule vibrational modes. Red, purple, blue and orange 
lines correspond to m = 0.1.3.5, respectively, where m represents the bosonic occupation of the 
vibrational modes. Standard paramenters are (0 = 800cm 1 , E\ — £ 2 — 1042cm ', A = 92cm 1 
and T 2 = 300K.^^ 


system can be enhanced by the exciton-vibration dynamics, as quantified by the trapping flux. In 
the regime of strong exciton-vibration coupling, namely, A 3> 1, the states become mainly vibra¬ 
tional and the excitonic transport is suppressed, as relfected by the decay of the flux illustrated in 
Figl2la), though the ETE still saturates to a perfect value of 100%. In this sense, we know that 
the energy transfer efficiency is necessary, but not sufficient for the description of the excitation 
energy transport in the photosynthesis antenna. In addition to this, other quantities, such as output 
work (heat current into RC) and EPR, are essential as well to completely describe the quality of 
this QHE on excitation energy transport. 

Moreover, by comparison of fluxes and ETE for different occupations of vibrations, it is clear 
that the higher excited mode of intra-molecule vibrations is more efficient for the enhancement of 
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ETE and flux trapped by RC. However, the environment around the chromophores always satisfies 
the condition Tico ~ k^Ti, thus only those lowest vibrational modes can be populated. Hence we 
show the results for m = 0,1,3,5 here. Particularly, the ETE and flux for m = 3 are of 72% and 
75% improvements with respect to m — 0, at A = 0.6. 

The energy transfer between the PEB 50 dimer and RC acquires the delocalization of the exci- 
tons. We now investigate the trends of delocalization under the influence of exciton-vibration inter¬ 
action. As is shown in Figj2^c), population of high-lying exciton state has a rapid, non-exponential 
decay as the coupling to vibrational modes increases, which can be traced back to the coherent 
transition from |l,0,m) to |0, l,m+ 1) and from |l,0,m+ 1) to |0,1. in). On the other hand, the 
delocalization of the wave packet is demonstrated in the large figure in Figj2jc), by increasing 
the coupling strength to vibrations. In fact, the coupling between excitons and vibrational modes 
leads to an effective amplification of electronic coupling by a factor of y/m+ 1 as mentioned be¬ 
fore, between the two molecules in a dimer. This is particularly reflected in the matrix elements: 
(l,0,m\H\0 ,1 ,m + 1), (1,0,m|//|0,l,m) and (1,0 ,m+ 11//|0,1 ,m). 


Coherence effect 

Since the dynamics of intra-molecule vibrations is considered together with excitons, we need 
to study the influence of exciton-vibration interaction on the coherence effect. Notice that the 
coherence effect mentioned here refers to the entanglement entropy and the sum of the off-diagonal 
elements of density matrix with different electronic states (quantum coherence). Due to the strong 
correltaion between entanglement and coherence, we use the quantum entanglement to quantify the 
coherence effect hereafter and the coherence will be enclosed in SI. To calculate the entanglement 
entropy, we first need to diagolize the density matrix at steady state p — Y,fPf\Vf){Wf\> where 
|t iff) = 'Lf%n\n\ (J') .niij') ■in(J')) and '?/ is the unitary transform matrix in the diagonalization of 
density matrix. For each component of pure state, the density matrix reads, by partial tracing over 
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the freedoms of other sites except the first one ( n 3 = m) 
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£(n 2 (v),m(v)|vA/)(^/|« 2 (v),m(v)) 

v=0 


7 / 7 3 

H I H FI (k) ,n s (j) Sn s ( k ),«,(/) 

y,/=0 \ ' k=0s=2 


niU))(M l )\ 


( 10 ) 


f (1) 

whose eigenvalues are A{_. Then the entanglement entropy of each pure component is Sj — 

+ hllnX'L), which subsequently gives the total entanglement entropy at steady state 


5 = 



(Afln^ + AW) 


( 11 ) 


Actually the formal definition of the entanglement entropy in mixture ensemble is still an open 
question. The rational for our choice of the expression in Eq.(??) for our system is: the exci- 
tonic basis describing the delocalized excitons is the one which gives the unique diagonal form 
of the steady-state density matrix and these delocalized excitons contain the coherence between 
the local sites and play an important role in experiments. Moreover, the off-diagonal elements of 
the Hamiltonian also leads to the entanglement in the excitonic basis. As illustrated by Fig|2td), 
the exciton-vibration interaction improves the coherence at first, owing to the energy-exchange 
between excitons and vibrational modes. However, based on the discussion on the delocaliza¬ 
tion above, we know that the interaction with vibrational modes leads to the enhancement of the 
delocalization of exciton wave packet. But we can further see from FigQc) that the population 
on site 1 of the excitonic state at large exciton-vibration coupling is closed to 0. This indicates 
the suppression of the back-and-forth oscillation between the two localized excitonic states, which 
consequently causes the downhill trends of the coherence. On the other hand, the exciton-vibration 
leads to the Rabi oscillation with the frequency ~ \/ 800 2 + A 2 ft) 2 , which indicates a rapid oscilla¬ 
tion as well as the mismatching between the energy splitting of excitons and quanta of vibrational 
modes at large coupling strength. Thus the quantum entanglement and coherence will eventually 
vanish in long time limit. This is partially in contrast to the previous predictions .- 17,31 Furthermore, 
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by the comparison between Figj2£a), |2£b) and[2ld), it is found that in the range of coupling to vi¬ 
brations before the extremum of entanglement, the improvement of quantum coherence gives rise 
to a rapid growth of flux trapped by RC and ETE, namely, A ~ 0.39,0.46,0.62,0.8 for m = 5,3,1,0 
respectively. In other words, the coherent energy transfer can lead to a significant improvement of 
quantum yield. The incoherent energy transfer becomes important on the further improvement of 
ETE after the coherent process. 


Macroscopic transport 

According to quantum thermodynamics, the observables on macroscopic level serve as an impor¬ 
tant role in the function of photosynthetic organisms as a QHE. We now investigate two represen¬ 
tative quantities: output work and entropy production rate (EPR). The former provides a critical 
measure of the macroscopic quality of the QHE and the latter quantifies the nonequilibriumness of 
the system on macroscopic level.The 1 st and 2 nd laws in thermodynamics give Q\ —Qi — Qtr — E, 
a + S — S t , where the output work to RC can be calculated by Q tr — Tr [H3> trap (p)\. S and S, are 
the rates of system entropy and EPR, respectively. Due to the assumption of large environments 
with the negligible back influence from system to environments, the entropy flux flowing from 
system to environments reads <7 = — + ^ + where T R c is the temperature of RC. At steady 

state, no energy consumption occurs inside the system so that Q tr — Q\ — Qi, S t = a. In several 
natural light-harvesting antennae the temperature in RC is always the same as that for low-energy 
solvated protein fluctuations, namely, around the room temperature. In Fig|2te) and[2]T), the out¬ 
put work and EPR have similar behaviors: the coupling of excitons to vibrational modes promotes 
the coherent energy transport which leads to a rapid growth of the work and EPR. The incoherent 
transport caused by higher coupling strength to vibrations further saturates the work and EPR to 
particular values (distinguished by different levels of intra-molecule vibrations). For the strong 
exciton-vibration coupling that A 1, the incoherent transport leads to the suppression of output 
work as well as EPR, since the vibration dominates the feature of states. 
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Effect of thermal relaxation on coherent energy transfer 


Now we study the influence of the low-energy thermal bath described by a continuous distribution 
of harmonic oscillators on the coherent energy transfer in the exciton-vibration dimer. The inter¬ 
action between exciton and bath is governed by Debye spectral density with the cut-off frequency 
(Dd — 50 fs -1 < ksT 2 . Here we consider two regimes of coupling strength between PEB 50 dimer 
and bath: very weak coupling Er = 4cm 1 and intermediate coupling Er = 34cm 1 . As is shown 
in Figj3] the dephasing and fluctuation induced by the low-energy protein motion gives rise to the 
promotion of quantum yield illustrated in Figj3]», since the motion of excitons globally becomes 
more delocalized when the coupling increases, as reflected in Fig j3fb). 




Coupling strength A to vibrational modes 


0.0 0.2 0.4 0.6 0.8 1.0 

Coupling strength A to vibrational modes 


Figure 3: (Color online) Comparison between the regimes of weak and intermediate couplings of 
exciton-bath, as a function of coupling to vibrational modes. Blue and purple curves correspond to 
the reorganization energies Er = 4cm 1 and 34cm -1 , respectively. 


Summary and remarks 

In conclusion, we investigated the effect of intra-molecule vibrational modes on the coherent en¬ 
ergy transfer in the light-harvesting complexes. It was demonstrated that the exciton-vibration 
interaction led to a non-trivial improvement of coherent energy transfer by the enhancement of 
coherence. Furthermore we also show that the exciton-vibration coupling can give rise to a perfect 
quantum yield (over 90%) for energy transport, at steady state. Our study provides the insights 
for the exploration of the intre- and intra-molecular vibrations on multi-molecule systems, i.e., 
the natural light-harvesting complexes FHCII with 32 chromophores and FMO complex with 7 
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molecules, to see how the exciton-vibration interaction affects the energy transfer pathways. 
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